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In  this  paper  an  algorithm  la  presented  which  solves  the 
linear  programming  problem.  This  algorithm  oombines  the  usual 


phase  one  (getting  feasibility)  and  phase  two  (getting 
optimality)  of  the  Simplex  or  Dual  methods  into  a  single 
phase.  The  algorithm  begins  with  either  a  single  activity 

(column)  or  a  constraint  (equation)  and  prooeeds  to  add 

*  »  » 

either  activities  or  constraints  one  at  a  time#  solving  the 
subproblems  which  arise  for  their  optimal  solutions.  The 
final  solution  is  attained  after  adding  the  last  activity  or 
constraint.  The  algorithm  promises  to  be  an  efficient  one  and 
has  several  advantages  which  arise  from  the  information  supplied 
about  subproblems. 
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SUItUlRY 

In  this  paper  an  algorithm  Is  presented  which  solves  the 
linear  programming  problem.  This  algorithm  combines  the  usual 
phase  one  (getting  feasibility)  and  phase  two  (getting 
optimality)  of  the  Simplex  or  Dual  methods  into  a  single 
phase.  The  algorithm  begins  with  either  a  single  activity 
(column)  or  a  constraint  (equation)  and  proceeds  to  add 
either  activities  or  constraints  one  at  a  time#  solving  the 
subproblems  which  arise  for  their  optimal  solutions.  The 
final  solution  is  attained  after  adding  the  last  activity  or 
constraint.  The  algorithm  promises  to  be  an  efficient  one  and 
has  several  advantages  which  arise  from  the  information  supplied 
about  subproblems. 
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THE  CROSS-SECTION  METHOD 
An  Algorithm  for  Linear  Programing 


INTRODUCTION 

Both  the  Simplex  [2]  and  Dual  [6]  methods  for  solving 
linear  programming  problems  require  two  distinct  phases.  In 
the  first  phase  a  feasible  solution  to  the  problem  is  found 
and  in  the  second  phase  an  optimal  solution  is  secured.  This 
paper  develops  a  method  for  solving  linear  programming  problems 
which  combines  these  two  phases.*  Por  this  and  other  reasons, 
which  are  examined  in  detail  in  Section  4,  it  promises  to  be 
more  efficient  than  existing  techniques.  Furthermore,  it  has 
certain  advantages  which,  apart  from  efficiency,  make  it  superior 
to  the  Simplex  or  Dual  methods  as  now  applied.  These  aa vantages 
are  due  to  the  fact  the  method  supplies,  in  addition  to  the 
final  answer,  the  optimal  solution  to  subproblems  constructed 
by  ignoring  certain  constraints  (equations)  or  activities 
(columns)  of  the  original  problem.  ThiB  characteristic  of  the 
technique  permits  one  to  compute  without  difficulty  solutions 
to  problems  which  differ  only  in  certain  activities  or  constraints, 
without  solving  separate  linear  programming  problems.  These 
advantages  are  described  in  Section  5* 

The  central  contribution  of  this  paper  is  a  method  of 
going  in  one  "step"  from  an  extreme  point  (e.p.)  feasible  dual 


# 

The  author  is  indebted  to  Dr.  Oeorge  B.  Dantzig  for 
posing  a  problem  which  inspired  this  investigation. 
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solution*  of  on*  linear  programing  problem  to  an  e.p.  feasible 
dual  solution  of  another  linear  programing  problem  consisting 
of  the  first  problem  with  an  additional  constraint.  This  leads 
to  a  technique  in  which  an  n-equation  problem  is  treated  by 
optimizing  in  successively  larger  (more  constrained)  subproblems. 
This  technique  is  described  in  Section  3*  The  Dual  Method  is 
used  to  achieve  optimization  in  each  subproblem.  Each  time  it 
is  achieved  for  a  subproblera,  another  constraint  is  brought  in 
by  our  method,  and  the  Dual  Method  is  again  applied.  The 
constraints  are  introduced  in  such  a  way  as  to  obtain  a  better 
bound  on  the  value  of  the  total  n-equation  problem  (like  an 
iteration  of  the  Dual  Method,  the  step  by  which  the  constraint 
is  introduced  drives  the  current  value  toward  the  optimal  one). 

In  the  first  section  we  develop  the  geometric  significance  of 
our  method  of  going  from  one  e.p.  feasible  dual  solution  to 
another  while  adding  constraints.  In  the  second  section  we 
derive  the  result  algebraically. 


1 .  0B0METHIC  MOTIVATION 

For  concreteness,  we  shall  always  deal  with  a  problem 
requiring  minimization.  Thus,  in  the  dual  problem  we  shall 


*  /  \ 

In  terms  of  (1)  below,  an  extreme  point  feasible  dual 

solution  may  be  considered  to  be  a  set  of  numbers  (t,,...,tt  ) 

such  that  for  m  indices,  i,,...,i  ,  we  have 

l  m 


•  * 


0 


k 


1,2,..., m. 


while  for  the  other  n  -  m  indices,  the  above  equation  holds  with 
"less  than  or  equal”  replacing  the  equality  sign. 
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always  maximize.  Our  method  requires  that  the  primal  problem 
be  in  equation  (as  opposed  to  inequality)  form.  This  is  no 
restriction,  since  "slacks”  can  always  be  added  to  a  problem 
to  put  it  in  that  form.  Hence  we  take  the  problem  to  be  that 
of  (1). 


allxl  +  a12x2  + 


(1) 


a21xl  +  a22x2  + 


amlxl  +  an2x2  + 


min  c^x^  +  c2x2  + 


+  amxn'bl  xi^0'  1 

+  a2nxn  *  b2 


+  *m*n  -  bm 


+ 


We  define  the  "kth  Bubproblem"  or  "k-problem,"  1  £  k  £  n, 
to  be  the  linear  programming  problem  which  arises  from  (l)  by 
removing  the  last  m-k  constraints.  We  call  an  e.p.  feasible 
dual  solution  to  the  kth  subproblem  a  ”kth  solution.” 

The  method  may  now  be  indicated.  Oiven  a  kth  solution, 
we  attempt  to  find  a  k+lst  solution  in  the  following  way. 

Since  the  kth  solution  corresponds  to  a  k— column  basis  (or 
kth  basis)  and  the  k+lst  solution  will  certainly  require  a 
k+l~column  basis*  we  attempt  to  find  a  single  column  which, 
added  to  the  columns  of  the  kth  solution,  will  provide  a 
k+lst  solution.  If  we  assume  for  concreteness  that  the  first 
k  columns  oa  (l)  were  in  the  basis  of  the  kth  solution,  we 
may  then  be  said  to  search  for  a  column  i,  i  >  k,  which  will 
make  the  basis,  indicated  in  (2)  correspond  to  a  k+lst 
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solution.  In  fact,  there  will  always  exist  such  an  index  i 
(assuming  there  tx'sts  a  solution  to  (l)).  Geometrically 
this  can  be  appreciated  as  follows!  Consider  a  2nd  solution 
to  (1),  and  let  m  -  4  in  that  problem.  This  2nd  solution 
corresponds  to  an  extreme  point  of  the  convex  polyhedron 
Indicated  in  Pig.  1.  The  two  lines  on  which  the  point  lies 
correspond  to  the  columns  of  the  2nd  basis  (for  definiteness 
we  take  these  to  be  the  first  and  second  columns)  and  the 
point  is  on  the  appropriate  side  (determined  by  the  inequality 

ifi  Tiau  l  cj 

and  indicated  by  an  arrow)  of  each  of  the  other  lines.  When 
the  third  equation  is  added  to  the  problem,  the  lines  become 
planes  and  in  particular  the  lines 

Tlall  +  T2a21  "  C1 
7rla12  +  r2a22  “  c2 
of  our  2nd  basis  become  the  planes 
TlaU  +  T2a21  +  ir3a31  ■  el 


Tla12  + 


2*22  +  t2*32  -  Cj  . 
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The  fact  that  wo  have  decided  to  leave  the  two  columns  of  the 
2nd  basis  in  the  3rd  basis  Indicates  that  our  3rd  solution 
will  lie  somewhere  on  the  intersection  of  the  two  planes  which 
the  first  two  columns  now  determine.  Our  second  solution  is 
clearly  of  this  type  but,  though  feasible,  it  is  no  longer  an 
extreme  point.  Hence  we  move  this  point  up  or  down  along  the 
intersection  of  the  two  planes  (we  discuss  later  whether  up 
or  down)  until  we  meet  a  third  plane.  (Our  use  of  the  word 
"move"  is  purely  rhetorical  since  the  change  will  not  require 
a  search  but  only  a  single  decision.)  This  determines  an 
extreme  point  and  becomes  our  3rd  solution. 


an  extreme  point 
feasible  dual  solu¬ 
tion  to  the  2-problem 


Pig.  1 


2.  ALOEBRAIC  DERIVATION 

We  assume  that  we  have  a  kth  basis  (corresponding  to  a 
kth  solution)  and  that  we  are  looking  for  an  index  i  which 
will  determine  a  k+lst  basis  of  the  form  assumed  in  (2).  The 
dual  variables  corresponding  to  the  kth  basis  are 
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(3)  ( •  •  •  ,vk)  "  *  *  *  ,ak^ 

and  by  assumption 


(c.  > c 


1*2' 


-  ej  i 0 


J  ■  1# 2 1 • • • »n* 


We  wish  to  choose  on  index  1  which  satisfies  the  condition 
analogous  to  (4)  for  the  k+l*>t  stage,  or 


(5) 


•#ck#ci) 


% 

-1 

/•u  \ 

all'  '  *  * ,alk'  ali 

•  •  * 

1  *2J  \ 

•  t  • 

•  1  • 

*kl'  *  *  * ,akk#  ^ 

• 

• 

i  'Vj  1 

*lc+l,l'  '  ••'ak+l,kak+l,l/ 

W+l.j/ 

J  -  1.2 

,n. 


By  straight forward  calculation  the  reader  may  verify  that  the 
inverted  matrix  involved  in  (5),  denoted  1b  related  to 

as  indicated  in  (6). 

'•  \ 

!  : 


z1#...,zk  t 


P-1490 

9-16-58 

-7- 


yl  /all 

:  ’  (  : 

yk  *ki 


(z1,...,zk)  -  fc(*k+l,l'“,J*k+l,k^^3' 


t  -  - 


a 


(Rk+l,l'***'ak+l,k^Bk1M  j  ~ ‘Sc+l,! 


li 


ki 


Applying  (6),  we  can  write  (5)  in  the  form 


)  ( ,  •  •  •  ,  Cj^i ) 


/ 


0  $  •  •  •  $  0 


°\ 
t 

0 


. 

yi\ 

• 

• 

• 

yk 

\z^,  ♦  •  • 

0/ 

\0, • • *|0 


°'l  ru 


*'J 


rcji°- 

% 


Distributing  the  matrix  multiplication  over  the  three  interior 
matrices  and  collecting  the  —  c^  with  the  first  term,  we  get 
the  following  throe  sunmands,  after  substituting  some  of  the 
values  given  in  (6)1 


(8) 


a 


(c1,...,ck)(Bk1) 


U 


®kj 


—  c 


J 


a 


li' 


—  ( c^, . . . , c^) (B^  )  .  J (z^, • • • $Z^) 

*ki 


1J 


a 


(9)  o^Zj, . .  .,zk)  .  J  -  j(clf  . .  •»ok)(Bk1) 

‘  -<a  ’  Sd 


11 


do) 


tak+l,jV 


Combining  (8),  (9)#  and  (10)  and  then  substituting  the  values 
for  (z  »^2*  *  *  *  ,zk)  and  *  given  in  (6),  we  have  this  equivalent 
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form  of  (7) i 


Thus  the  oondition  to  be  satisfied  is  of  the  form 


(12) 


cj  1  ^ 
J  a 


- 

—jL  c 


J  ™  1|2# • ■ • 


where  the  are  relative  cost  factors  for  the  kth  subproblera 
and  the  a,,  J  -  l,...,n,  are  the  numbers  which  would  be  the 
relative  cost  factors  for  the  kth  subproblem  if  the  coefficients 
of  the  next  equation  were  taken  as  the  cost  form  for  the  problem. 
Since  B^.  is  a  kth-basis,  the  c y  J  «  l,...,n,  are  all  non¬ 
positive.  Using  this  fact,  the  reader  can  verify  that  if  there 
exists  an  >0,  then  an  index  i^  for  whioh 

c.  Ci1 

a.  >0  and  — £  — —  for  all  J  such  that  a«  >  0 
1  8i, 


(13) 
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wlll  satisfy  (12).  Similarly,  If  there  exists  an  a^  <  0,  then 
an  Index  i2  for  which 

C .  cl„ 

(14)  a.  <0  and  -*•  >  -  for  all  J  such  that  a.  <  0 

2  %  J 

will  satisfy  (12).  We  are  precluded  from  choosing  an  Index  1 
for  which  -  0  since  (12)  would  not  be  defined.  This 

corresponds  to  the  fact  that  the  addition  of  such  a  column 
makes  a  singular  matrix  (without  the  necessary  inverse). 

With  this  understanding,  it  is  easily  verified  that  any  index  1 
for  which  (},2)  holds  must  satisfy  either  (13)  or  (l4).  We  now 
consider  the  implications  of  a  choice  between  (13)  and  (14). 

If  In  (5)  we  delete  c^  and  let  the  column  of  b's  in  (l) 
play  the  role  of  the  Jth  column  in  (5)*  the  left-hand  side  of 
that  equation  will  represent  the  value  associated  with  the 
k+lst  solution  whioh  we  denote  Reasoning  similarly  with 

(4)  we  can  interpret  the  reduction  of  (5)  to  ( 11 )  as  indicating 
that 


where  F  is  defined  analogously  to  a^.  Since  we  are  interested 
in  maximizing  in  the  dual  problem,  we  evidently  should 
ascertain  the  sign  of  F  and  choose  our  index  i  in  such  a  way 
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D  —  .  * 

that  —  is  positive.  Assuming  that  b  f  0,  this  is  always 

5i 

possible,  provided  that  the  equations  have  a  feasible  primal 
solution. 

The  reasoning  is  as  follows:  Let  us  pretend  that  the 
k+lst  equation  represents  the  cost  form  for  the  problem; 

F  becomes  the  difference  between  the  present  value  and  the 
value  which  must  be  achieved  to  satisfy  the  k+lst  equation. 

The  a^  become  the  relative  cost  factors.  If  E  (  0,  and 
2  0*  J  *  l,...,n,  we  would  have  a  feasible  dual  solution 
to  a  problem  of  maximization  (in  the  primal)  in  which  more 
"value”  was  required  to  satisfy  equation  k+1.  Since  the 
value  associated  with  the  feasible  dual  solution  is  an  upper 
bound  on  the  value  which  can  be  attained  by  a  feasible  primal 
solution,  there  can  exist  no  feasible  primal  solution  to  the 
k+l-problem  (and  hence  no  such  solution  to  the  complete  m- 
problem) .  A  similar  argumer'  applies  to  the  case  in  which 
F  >  0  and  a<  £  0,  J  -  l,...,n.  This  gives  the  result  and 
allows  us  to  Bpy  that  the  bound  on  the  final  optimal  solution 
attained  by  bringing  in  the  equation  is  no  worse  (and  usually 
better)  than  the  previous  one.** 

#  __ 

If  d  «  0,  the  k+lst  equation  is  satisfied,  along  with  the 
first  k  equations,  by  variables  which  correspond  to  a  k-column 
basis.  In  this  case  any  column  may  be  brought  in  which  satisfies 
a^  /  0.  If  no  such  oolumn  exists,  the  k+lst  equation  is  a  linear 
combination  of  the  previous  k  equations  and  can  thus  be  ignored. 

*  •» 

Dr.  Oeorge  Dantzig  has  pointed  out  that  the  result  in 
this  section  bears  a  close  relationship  to  the  Parametric 
Linear  Programming  of  William  Orchard-Hays  [7]  .  Looked  at 
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3*  ALOORITHM  PRBSSNTATION 

Our  method  of  Introducing  equations  suggests  an  algorithm, 
which  will  be  discussed  in  this  section,  for  solving  linear 
programming  problems.  We  propose  to  solve  the  »-equation 
problem  by  the  following  steps i  Beginning  with  an  (e.p.) 
optimal  solution  to  a  one— equation  eubproblem,  we  will  introduce 
another  constraint.  This  can  be  done,  because  an  e.p.  optimal 
solution  to  a  problem  corresponds  to  an  e.p.  feasible  dual 
solution.  Our  method  derives  an  e.p.  feasible  dual  solution 
to  the  2-problera,  and  using  the  Dual  Method  we  aohieve  an 
optimal  solution  to  this  subproblem.  We  are  then  in  a  position 
to  introduce  another  constraint,  giving  a  three-equation 
subproblem.  We  continue  in  this  way  until  we  have  brought  in 
the  last  equation  and  used  the  Dual  Method  to  find  an  optimal 


(cont'd  from  page  10) 

in  this  light,  the  application  of  the  method  corresponds  to 
moving  the  right-hand  side  of  a  k+lst  subproblem  "parametrically” 
from 


Our  approach  differs  from  Parametric  Linear  Programming  in  that 
we  do  not  necessarily  begin  with  an  optimal  situation  and  do 
not  seek  to  maintain  optimality.  Further,  the  movement  indicated 
above  is  carried  out  not  "parametrically”  but  in  one  step. 
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solution  to  that  ultimata  subproblem. 

There  is  only  one  difficulty  to  to  considered  with  respect 
to  this  procedure.  It  is  possible  for  all  of  the  m  possible 
one-equation  subproblems  to  fail  to  have  an  (e.p.)  optimal 
solution  although  the  entire  a-equation  problem  is  quite 
unexceptional  (i.e.,  has  optimal  programs).  In  this  case  we 
employ  a  trick  due  to  Beale  [l]  and  begin  by  considering  a 
fictitious  equation, 

(16)  Z  x.  +  t  -  M, 

i-1  1 

where  M  is  a  very  large  positive  number  (thought  of  as  larger 
than  the  sum  of  the  x^  for  any  feasible  primal  solution,  though 
in  fact  there  may  be  no  such  nunber) .  If  we  assign  t  a  zero 
co jt  this  1-problem  has  an  (e.p.)  optimal  primal  solution. 


# 

We  could  have  proposed  In  place  of  this  procedure  that 
the  dual  of  the  original  problem  be  considered  the  primal 
problem  and  the  above  procedure  followed.  This  would  lead  to 
maintaining  feasibility  in  the  original  problem  and  successively 
adding  activities  (columns).  This  important  possibility  is 
discussed  in  Section  5. 

•  • 

An  example  is  shown  below.  Each  equation  in  devoid  of 
feasible  dual  solutions  as  a  1-problem  but  together  they  are 
not.  An  optimal  solution  is  x^  -  0,  xn  -  2,  -  1. 

X1  *2  x3 

r1  1  -1  0  -  -2 

1  -1  -1  - 

rain  -1  2  -1 


-3 
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and  hence  we  are  in  a  petition  to  initiate  our  procedure.  The 
addition  of  (16)  to  the  original  problem  does  not  diminish 
(hence  does  not  change)  the  set  of  feasible  primal  solutions 
to  the  problem  if  this  set  is  bounded.  It  does,  however# 
ensure  that  the  set  of  feasible  dual  solutions  is  nonempty# 
and  hence,  by  the  Duality  Theorem  [41 ,  it  guarantees  the 
existence  of  a  solution  (which  our  algorithm  will  find)  to  the 
modified  problem.  If  the  dual  variable  associated  with  (16) 
ever  becomes  zero  (as  it  usually  will  early  in  the  procedure)# 
(16)  can  be  ignored  thereafter.  However#  it  is  possible  that 
in  the  optimal  solution  this  dual  variable  is  not  zero.  This 
indicates  that  the  set  of  values  of  feasible  primal  solutions 
to  the  original  problem  is  not  bounded  below  (i.e.#  the 
original  problem  haul  an  infinite  solution) . 

4.  EFFICIENCY  IN  COMPUTATION 

We  turn  now  to  considering  the  efficiency  of  our  proposed 
algorithm.  If  we  were  to  introduce  one  equation  after  another 
without  attempting  to  find  optimality  in  each  subproblem#  we 
would  eventually  reach  an  e.p.  feasible  dual  solution  to  the 
entire  problem.  Thus  our  result  would  be  similar  to  that 
attained  by  a  phase  one  algorithm  for  the  Dual  Method.  However, 
because  we  would  have  introduced  constraints  in  such  a  way  as 
to  increase  the  value  of  the  dual  solution,  we  might  expect  to 

If  the  artificial  equation  is  required,  the  e.p.  feasible 
dual  solution  will  be  associated  with  a  slightly  more  constrained 
problem. 
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be  somewhat  nearer  the  optimal  solution  than  would  be  the  case 
if  we  had  followed  a  phase  one  procedure.*  Neglecting  this 
advantage,  we  might  consider  the  introduction  of  equations 
in  our  algorithm  to  be  analogous  to  a  phase  one  technique. 

Hence  we  might  compare  the  number  of  multiplications  required 
to  introduce  the  entire  set  of  equations  by  our  method  with 
the  number  of  multiplications  required  by  some  phase  one 
procedure.  The  number  of  multiplications  we  require  is,  of 
course,  not  dependent  upon  whether  we  pause  to  optimize  in 
the  subproblems  or  Just  pass  quickly  from  one  to  the  next. 

Having  made  the  above  comparison,  we  will  then  go  on  to 
compare  the  number  of  multiplications  which  might  be  necessary 
to  reach  optimality  in  all  the  subproblems  with  the  number 
usually  required  to  perform  a  phase  two  (i.e.,  reach  optimality, 
assuming  feasibility)  with  the  Dual  or  Simplex  Method. 

At  the  k+lst  stage,  two  sets  of  k  dual  variables  must  be 
computed  in  our  algorithm,  one  set  corresponding  to  the  true 

cost  form  and  one  set  corresponding  to  the  pretense  that  the 

2 

next  equation  is  the  cost  form.  This  requires  2k  multiplications 


» 

Only  two  approaches  to  the  problem  of  constructing  a 
phase  one  procedure  for  the  Dual  Method  are  known  to  us. 

Lemke  in  [b]  proposes  that  a  Simplex  Method  phase  two  be 
performed  on  a  problem  which  he  describes.  0.  Dantzig  in  [3, 
suggests.  Instead  of  a  phase  one  for  the  Dual  Method,  that  the 
Dual  Method  be  applied  to  a  problem  related  to  the  original 
one  for  which  an  e.p.  feasible  dual  solution  is  insnediate  and 
that  the  Simplex  Method  be  used,  after  the  Dual  Method  has 
found  the  optimal  solution,  to  go  from  the  optimal  solution 
of  the  related  problem  to  the  optimal  solution  of  the  original 
problem. 
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in  all.  Applying  these  sets  of  dual  variables  to  the  (n-k) 
columns  not  yet  in  the  basis,  in  order  to  compute  the  and 
a<,  requires  2k (n-k)  multiplications.  Determination  of  F 
requires  k  more  multiplications,  and  a  number  of  ratios  not 
exceeding  n-k  must  then  be  computed. #  Hence  the  transition 
requires  less  than 

(17)  2k2  +  2k(n-k)  +  k  +  n  -  k  -  2nk  +  n 

multiplications.  Sumning  (17)  as  k  goes  from  0  to  »-l  will 
indicate  the  number  of  multiplications  necessary  to  bring  in 
all  a  equations.  Hence  we  have 
m-1 

(18)  2  2nk  +  n  -  2n 
k-0 

Counting  nm  multiplications  for  a  Simplex  or  Dual  Method 
Iteration,  we  interpret  this  result  as  the  equivalent  of 
m  iterations.  It  seems  reasonable  that  no  phase  one  procedure 
can  be  consistently  faster  than  this.  For  example,  it  would 
require  at  least  m  iterations  to  remove  m  artificial  vectors, 
if  these  were  used  in  a  phase  one  procedure.  A  phase  one 
procedure  that  Involved  an  auxiliary  Simplex  or  Dual  problem 
would  also  require  at  least  m  iterations  most  of  the  time. 


P-1490 

9-16-58 

-l6- 


If  we  nay  Judge  by  computing  experience.  We  conclude, 
therefore,  that  the  Introduction  of  equation*  by  our  algorithm 
compares  favorably  in  speed  with  a  phase  one  technique  as 
well  as  having  the  advantage  mentioned  at  the  start  of  this 
discussion. 

The  comparison  of  a  phase  two  technique  with  the  Dual 
Method  iterations  necessary  for  our  algorithm  to  reach 
optimality  in  every  subproblem  is  best  carried  out  in  geometric 
language.  Geometrically  our  proposal  reduces  to  considering 
certain  cross— sections  of  the  polyhedron  of  feasible  dual 
solutions.  At  the  kth  stage  the  cross-section  considered 
lies  in  the  coordinate  "face"  determined  by  «  0, 
i  -  k+1,  k+2,  ...»  m.  (In  other  words,  the  cross-section 
considered  is  the  one  achived  by  ignoring  the  last  m-k  equations.) 
Our  algorithm  then  uses  the  Dual  Method  to  optimize  in  that 
cross-section.  When  no  better  bound  on  the  optimal  value  can 
be  achieved  in  the  cross-section,  another  equation  is 
introduced  (this  operation  also  improves  the  bound)  and  the 
Dual  Method  is  again  applied.  The  cross-section  now  lies  in  a 
coordinate  face  determined  by  -  0,  I  -  k+2,  k+3,  m. 

In  this  "higher  dimensional"  crosa-aecticn  our  algorithm 
prohibits  us  from  considering  any  "vertices"  (of  the  cross- 
section)  whioh  are  associated  with  lower  values  of  the  dual 

This  statement  is  appropriate  as  a  comment  on  Leake's 
method  for  constructing  a  phase  one  procedure  isee  previous 
footnote).  It  is  difficult  to  determine  the  effect  of  0.  B. 
Dantzig's  proposal. 
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than  have  been  already  achieved.  Hence  we  consider  a  truncated 
cross— section.  It  is  our  expectation  that  achieving  optimality 
for  the  previous  subproblem  and  introducing  the  next  equation 
as  we  do  (so  as  to  achieve  a  better  bound)  will  truncate  the 
next  cross-section  radically,  so  that  few  iterations  will  be 
required  to  regain  optimality.  This  expectation  is  one-half 
of  the  argument  that  our  algorithm  provides  a  good  alternative 
to  a  phase  one  technique.  The  other  half  of  our  argument 
centers  on  the  fact  that,  whatever  the  number  of  iterations 
required  to  achieve  optimality  for  eaoh  subproblem  may  be, 
these  iterations  often  concern  a  relatively  email  basis  and 
hence  involve  fewer  multiplications  than  a  Simplex  or  Dual 
Method  iteration.  For  example,  let  us  consider  the  possibility 
that  it  requires  two  Dual  Method  iterations  at  eaoh  stage  to 
achieve  optimality  for  the  subproblem.  In  this  case,  counting 
kn  multiplications  for  a  Dual  Method  iteration  in  a  k-equation 
problem  we  would  require 

Z  2kn-  gnm(m  + 1)  _  2n  «  nm(m  +  1)  -  2n 
k-2  '  * 

multiplications.  This  result  is  somewhat  less  than  the 
equivalent  of  m+1  Simplex  or  Dual  Method  iterations.  With 
the  m  iterations  required  to  introduce  the  equations  included, 
our  entire  algorithm  would  take  the  equivalent  of  approximately 
2m  Simplex  iterations  and  would  probably  be  about  on  a  par 
computationally  with  the  techniques  in  use. 

It  3b  obviously  possible  that  it  will  take  somewhat  more 


than  two  Dual  Method  iteration*  to  achieve  optimality  for  each 
subproblem.  It  is  relevant  to  point  out  that  the  new  variable 
x^  introduced  into  the  basic  solution  at  the  k+lth  stage 
always  enters  at  a  nonnegative  level.  It  is  derived  in  a 
straightforward  manner  to  be 


(19) 


a 


i 


which  is  nonnegative  since  is  always  chosen  to  be  of  the 
same  sign  as  t>.  Of  course  the  k+lst  solution  may  not  be 
optimal  when  the  kth  solution  is,  because  any  of  the  first  k 
variables  may  become  negative.  The  now  values  of  the  first  k 
variables  (primed)  have  the  relationship  to  the  old  values 
indicated  in  (20). 


(20) 


; 


In  any  case,  the  nonnegativity  of  x4  allows  for  the  possi¬ 
bility  that  as  the  equations  are  brought  in  an  optimal  solution 

is  maintained  from  the  kth  to  the  k+lst  stage  (without  the 

\  * 

use  of  any  Dual  Method  iterations). 


A  related  question  to  which  the  answer  is  definitely  in 
the  negative  askB  whether  there  always  exists  some  ordering  of 
the  equations  which  will  allow  them  to  be  brought  in  one  by 
one  while  the  kth  solutions  remain  optimal  at  each  stage.  A 
counter-example  is  given  below.  Depending  upon  whioh  equation 
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5.  ADVANTAGES  IN  APPLICATION 

In  this  section  we  discuss  the  flexibility  and  advantages 
of  the  algorithm  proposed.  We  begin  by  presenting  three  types 
of  situations  in  which  the  algorithm  would  prove  useful. 

A.  A  firm  has  listed  the  constraints  whioh,  it  feels, 
might  apply  in  a  given  problem.  It  recognizes  that  if  certain 
constraints  prove  too  restraininc  in  terms  of  the  optimal 
value,  they  might  be  altered  or  removed  by  some  "nonlinear  action, 
such  as  acquisition  of  more  land  or  plant  facilities,  etc.  Our 
recommendation  is  to  place  these  alterable  constraints  last,  and 
to  use  the  cross-section  method  to  achieve  an  optimal  solution 
to  the  subproblem  P,  consisting  of  the  other  k  constraints,  and 
thereafter,  for  each  alterable  equation,  to  alternate  bringing 
in  the  equation  by  our  method  and  reaching  optimality  in  that 
k+  1-equation  subproblem  by  the  Dual  Method,  in  order  to  gauge 
the  difference  between  the  optimal  value  of  p  and  the  optiral 
value  of  the  subproblem  consisting  of  P  with  an  alterable 

(cont'd  from  page  18) 

is  considered  initially,  column  one  or  column  4  is  introduced 
into  the  basis.  However,  neither  of  these  columns  appears  in 
the  only  optimal  solution  to  the  problem.  That  solution  is 

X1  *  0#  *2  "  V3#  x3  “  V3#  *4  -  0. 


*1 

*2 

x3 

x4 

T1  1 

1/2 

1/* 

1/8 

T2  1/3 

i/* 

1/2 

1 

min  —1 

-1 

-1 

-1 

1 
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equation  added.  This  will  indioate  the  restraining  effect  of 
any  single  alterable  equation  on  the  set  of  unalterable 
equations . 

B.  A  firm  has  listed  the  inequalities  which  it  feels 
might  apply  in  a  given  problem  but  recognizes  that  most  of 
them  are  probably  not  constraining  in  an  optimal  solution. 

Rather  than  work  with  the  entire  set  of  Inequalities,  as  would 
be  necessary  by  a  straightforward  application  of  the  Simplex 
Method,  the  firm  should  place  at  the  end  those  inequalities 
which,  in  its  Judgment,  are  thought  to  be  unimportant  and 
then  proceed  to  apply  the  technique  indicated  in  Section  3 
(after  adding  slack  vectors).  A  check  on  the  rate  at  which 
the  optimal  values  to  successive  subproblems  decrease  should 
provide  a  good  indication  of  how  restraining  the  latest  equations 
are.  When  the  values  cease  to  decrease,  the  firm  might  want 

to  assume  that  the  remaining  equations  are  not  constraining. 

In  this  case  it  would  have  only  to  check  the  feasibility  (in 
the  other  inequalities)  of  the  optimal  solution  to  the  last 
subproblem  considered.  Assuming  that  the  constraints  could 
be  satisfied,  the  firm  would  have  an  optimal  solution,  without 
the  labor  involved  in  handling  an  oversized  matrix. 

C.  A  firm  might  be  interested  in  a  bound  on  the  optimal 
value  of  a  given  linear  programming  problem.  It  might  feel 
that  if  the  optimal  value  v  were  greater  than  some  Vq,  it 
would  be  pointless  to  follow  some  course  of  action  associated 
with  the  problem.  In  this  case,  the  procedure  advised  in 
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Section  3  would  give  successively  better  lower  bound#  as  each 
8ubproblera  was  considered.  If  the  lower  bound  indicated  by 
some  subproblem  surpassed  Vq,  the  firm  would  know  that  the 
course  of  action  considered  could  not  be  followed  unless  at 
least  one  of  the  constraints  involved  in  that  optimal  sub- 
problem  were  altered  or  removed. 

The  advantages  discussed  so  far  exhibit  the  virtues  of  a 
method  which  adds  equations  one  at  a  time  while  completely 
solving  the  Induced  subproblems .  These  advantages  would  have 
their  counterpart  in  a  method  which  adds  activities  (columns) 
one  at  a  time  while  completely  solving  the  subproblems  which 
arise.  In  order  to  get  such  a  method  we  have  only  tc  apply 
the  procedure  we  recoimnended  in  Section  3  to  the  uual  of  the 
original  problem.  Suoh  a  procedure  maintains  feasibility  of 
the  intermediate  solutions  to  the  original  problem  and  henee 
provides  upper  boundB  (according  to  our  conventions),  as 
opposed  to  the  lower  bounds  utilized  in  example  C.  Furthermore 
it  allows  one  to  alter  activities,  with  all  the  advantages 
analogous  to  those  described  in  examples  A  and  B  for  ohanglng 
equations. 
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